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ABSTRACT 

The hydrogen ionization and dissociation front around an ultraviolet radiation source should merge when the 
ratio of ionizing photon flux to gas density is sufficiently low and the spectrum is sufficiently hard. This regime 
is particularly relevant to the molecular knots that are commonly found in evolved planetary nebulae, such as 
the Helix Nebula, where traditional models of photodissociation regions have proved unable to explain the high 
observed luminosity in H2 lines. In this paper we present results for the structure and steady-state dynamics 
of such advection-dominated merged fronts, calculated using the Cloudy plasma/molecular physics code. We 
find that the principal destruction processes for H2 are photoionization by extreme ultraviolet radiation and 
charge exchange reactions with protons, both of which form H2^, which rapidly combines with free electrons 
to undergo dissociative recombination. Advection moves the dissociation front to lower column densities than 
in the static case, which vastly increases the heating in the partially molecular gas due to photoionization of 
He", H2, and H°. This causes a significant fraction of the incident bolometric flux to be re-radiated as thermally 
excited infrared H2 lines, with the lower excitation pure rotational lines arising in 1000 K gas and higher 
excitation H2 lines arising in 2000 K gas, as is required to explain the H2 spectrum of the Helix cometary 
knots. 

Subject headings: hydrodynamics — molecular processes — planetary nebulae: individual (NGC 7293) 



L INTRODUCTION 

The ultraviolet photons from hot stars, such as main- 
sequence OB stars or the central stars of planetary nebula 
(CSPN) will dissociate and ionize the surrounding ckcumstel- 
lar and interstellar ga s. In the classical picture (e.g., |Hollen-| 
bach & Tielens|1997 1, the extreme ultraviolet (EUV) photons 
with energies hv > 13.6 eV photoionize the hydrogen gas, 
forming an H II region bounded by a relatively sharp ioniza- 
tion front (IF), while the far ultraviolet (FUV) photons with 
energies 6 eV < /iv < 13.6 eV penetrate the IF to form a neu- 
tral photodissociation region (FDR) that surrounds the H II 
region. The dissociation of H2 in such a PDR is dominated 



by a two-step radiative process (Stecher & Williams 1967 



Abgrall et al. 2000), in which absorption of an FUV photon 
leaves the H2 molecule in an electronically excited state, from 
which it has a certain probability 10%) of decaying to the 
vibrational continuum of the ground electronic state. 

However, as shown by Bertoldi & Draine ( 1996), a classi- 
cal extended PDR cannot exist if the FUV flux is sufficiently 
weak compared with the EUV flux at the IF, rather the IF and 
dissociation front (DF) should merge. Bertoldi & Draine con- 



sidered the case of H II regions around OB stars and found 
that this regime was most relevant for cases in which the dust 
optical depth through the ionized gas is low, which corre- 
sponds to a low ionization parameter (the ionization parame- 
ter is a dimensionless number equal to the ratio of the number 
density of ionizing photons to the number density of hydro- 
gen nuclei). To date, no models have been calculated of the 
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structure and emission properties of such fronts, which are 
also expected to show strong deviations from static chemical 
and ionization equilibrium. 

In this paper, we calculate in detail the internal dynam- 
ics and chemistry of merged ionization/dissociation fronts 
(IF/DFs), concentrating on the particular case of compact 
photoevaporating molecular knots in evolved planetary neb- 
ulae (PNe), such as are seen in the Helix nebula (|Young et aL 



1999 



IMeixner et aL]|2005l |Hora et al.||2006l [O'Dell et al. 



2007 OHF07). The stellar spectrum from the hot central star 
of a PN is harder than that from an OB star, leading to an 
EUV luminosity that is much greater than the FUV luminos- 
ity. In addition, the ionization parameter of the knots is much 
lower than is typically seen in H II regions, resulting in very 
little attenuation of the EUV flux by recombinations in the 
ionized gas. Both these factors place the knots firmly in the 
merged IF/DF regime. The most comprehensive existing the- 
oretical study of PDRs in planetary nebulae is that of Natta] 



& Hollenbach ( 1998) , who present detailed modeling of the 



time-dependent evolution of an expanding circumstellar shell 
as the luminosity and spectrum of the CSPN evolves, finding 
that soft X-rays can be important in exciting the molecular 
gas at late times. However, OHF07| showed that this is not 
the case in the Helix since it is only in the EUV band that the 
CSPN luminosity is sufficient to excite the knot PDRs. Natta] 
& Hollenbach used an analytic treatment of the H II region, 
which is assumed to absorb all the EUV radiation, and were 
thus unable to model the case of a merged IF/DF. 

2. MODELS 

In order to investigate the structure of advective IF/DFs 
in PNe, we have calculated steady-state, plane-parallel slab 



models using the Cloudy plasma/molecular physics code (Fer- 
land et aL||I998 1. Details of the computational scheme used 
to treat the steady-state dynamics are given in |Henney et al. 



(2005) and these metho ds have now been coupled to a hydro- 
gen chemical network ( Abel et aLpOOS^ and combi ned with 
the 1893-level model of H2 described in fshaw etaL] ( |2005] l. 
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TABLE 1 
Model input parameters 



Stellar parameters 










L = 120 Lo 


Teff = 


130,000 K 








Qh = 7.57 X 10"^ s- 


' Op, IV = 


1.33 X 10"*^ s"' 








D 


"0 


"0 


Fo 




Model 


(pc) 


(kms-') 


(cm-3) 


(cm^^ s"' ) 




AlO 


0.137 


10 


3162. 


3.36 X 10' 


0.94 


A06 


0.137 


6 


3162. 


3.36 X 10' 


0.56 


A03 


0.137 


3 


3162. 


3.36 X 10' 


0.28 


AOl 


0.137 
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3162. 


3.36 X 10' 


0.09 


AOO 


0.137 





3162. 


3.36 X 10' 


0.00 


AAIO 


0.137 


10 


1000. 


3.36 X lO' 


0.30 


BIO 


0.244 


10 


1000. 


1.06 X 10' 


0.94 


B06 


0.244 


6 


1000. 


1.06 X 10' 


0.57 


BOO 


0.244 





1000. 


1.06 X 10' 


0.00 


CIO 


0.433 


10 


316. 


3.36 X lO*' 


0.94 


C06 


0.433 


6 


316. 


3.36 X 10** 


0.56 


COO 


0.433 





316. 


3.36 X 10** 


0.00 



3. PREDICTED MODEL STRUCTURE 

Figure [T] shows the results from a typical advective model, 
B06, while Figure |2] shows results from a static model, BOO, 
with exactly the same incident flux and density as B06. For 
ease of discussion, we divide the model into 3 broad zones: 
Zone I is closest to the ionizing source and is largely ionized, 
with a very low molecular fraction; Zone II is the dissociation 
front, in which the hydrogen is half neutral and half molecu- 
lar (for the advective models, this zone is subdivided into Ila 
and lib); Zone III is farthest from the ionizing source, where 
hydrogen is fully molecular Table [2] shows typical physical 

* Note that all flow velocities, u, are in the frame of reference of the IF/DF, 
but since we find that u becomes very small at great depths, this is also ap- 
proximately the rest frame of the molecular gas. 
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Fig. 1. — Structure of a representative advective model (B06, see Table^ 
as a function of column density of hydrogen nuclei. Top panel: Various flow 
variables on a logarithmic scale: flow velocity, u, depth into cloud, z, to- 
tal number density of hydrogen nuclei, n, hydrogen molecular fraction, /h, , 
electron fraction n^jn, and gas temperature, T. Bottom panel: Emissivity of 
important coolant lines. 



In this initial study, we restrict oursel ves to models with 
elemental abundances from|Henry et al.|(]1999) that are illu- 
minated by a black-body source of luminosity 120 L© ™d ef- 
fective temperature 130,000 K, chosen to match the CSPN of 
the Helix nebula ( |BohUn et al.|1982] adjusted for the trigono- 
metric parallax of 217 pc, Harris et al. 20071. We have also 
calculated some models using a |Rauch ( 2003| stellar atmo- 
sphere with the same luminosity and efifective temperature. 

We vary three different model parameters so as to span the 
range of physical and illumination conditions that are seen 
in PN knots: distance from the CSPN, D; hydrogen nuclei 
density at the illuminated face, «o; and gas velocity at the il- 
luminated face, M(). Table [T] summarises the input parameters 
of our models.^ 

The magnitude of advective eff'ects in the models is, to 
first order, dependent only on the dimensionless combination: 
'tad - nouoAnD^ / Qu, where Qu is the ionizing photon lumi- 
nosity of the source. This advection parameter ( |Henney et al. 
[2005) is the ratio of particle flux to photon flux, which in- 
creases with Mo and decreases with ionization parameter. The 
appropriate value of the downstream flow velocity mq depends 
sensitively on the boundary conditions at the illuminated face 
and on the global geometry of the flow. In the case of a 
photoevaporating knot with negligible confining pressure on 
the ionized side, uq will be of order the ionized sound speed 
10 km s"'), and this is the case we concentrate on in this 
paper. In the case of a more shell-like configuration of the 
molecular gas, mq would tend to be lower. 
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Fig. 2. — Same Fig.^ but for an equivalent static model (BOO). 
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Fig. 3. — Destruction and formation rates of Ht, in units of molecules per 
hydrogen nucleus per second, for the advective model shown in Fig. [T] Top 
panel: Formation rates. Bottom panel: Destruction rates. Note the dnferent 
scales of the y axes. 

conditions in each zone for the advective models. 

The diff'erences between the advective and static models are 
stark: in the advective model, the DF occurs at the very low 
column density of 10'^ cm"^ from the illuminated face and 
strongly overlaps with the ionization front, whereas in the 
static model the DF occurs much deeper, at 2 x 10''^ cm"^. 



Merged ionization/dissociation fronts 
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TABLE 2 

Physical conditions in different zones of a typical advective IF/DF structure 



Zone 


Column (cm ^) 


r(K) 


u (km s ' ) 




/h. 


n/no 


Heat 


Cool 
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< 10'** 


10^ 
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w-^ 
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H" p.e. 


Metal, H" lines 
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Fig. 4. — Thermal balance of the advective model shown in Fig. [T] Top 
panel: Principal contributions to the atomic/molecular heating rate per unit 
mass, shown on a logarithmic scale. Middle panel: Same as top panel, but 
for cooling rate. Bottom panel: Difference between heating and cooling rates, 
shown on a linear scale. 

in a region where the ionization fraction is < 0.01. Zone II 
is much warmer in the advective model (2000 K in Zone Ila, 
1000 K in Zone lib) than in the static model (500 K) and, as 
a result, H2 line emission is more than an order of magnitude 
brighter 

The rates of formation and destruction of H2 for Model B06 
are shown in Figure|3] where it can be seen that the destruction 
rate (bottom panel) has a narrow peak at the leading edge of 
Zone Ila, due to collisional processes with protons and elec- 
trons, together with a broader peak covering Zones Ila and 
lib, due to photoionization by hard EUV photons. The prin- 
cipal reaction channel for the ions formed by these pro- 
cesses is dissociative rec ombination with free electrons (e.g., 
McCandhss et al.||2007| ), with only ~ 10% reacting with H° 
to re-form H2. Other H2 formation mechanisms have even 
smaller rates (top panel), with the result that, once they have 
been destroyed, most H2 molecules never re-form during the 
^ 300 yr that they remain in the DF. 

The heating and cooling rates for Model B06 are shown in 
Figure |4] In Zone I, as is typical for low -excitation H II re- 
gions, the heating is dominated by photoelectric heating, 
while the cooling is due to H lines and optical metal lines 
such as [N II] 6584 A. In Zone Ila, H° photoelectric heating 
still dominates the heating, whereas in Zone lib, photoelec- 
tric heating of He" and H2 increasingly take over In all of 
Zone II the cooling is overwhelmingly dominated by H2 line 
emission. In Zone III, the heating rate is much lower than in 
the other zones and is due principally to the absorption of H2 
lines emitted in Zone II, while the cooling in Zone III is dom- 
inated by coUisionally excited far-infrared metal lines. The 
bottom panel of Figure |4] shows the difference between the 
heating and cooling rates, which is equal to the net rate of en- 
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Fig. 5. — Fraction of the stellar bolometric luminosity that is radiated in 
H2 lines, assuming 100% covering fraction of molecular gas. Black symbols 
show the total H2 line luminosity, white symbols show the sum of the pure 
rotational 0-0 mid-infrared lines from 5(0) to S(28), gray symbols show the 
near-infrared 1-0 S(l) line at 2.121 pm. Different symbol shapes are for 
different sequences of models, as shown in the key. The dashed curve shows 
an approximate analytic fit (see text). 

ergy transfer from the radiation field to the gas. This can be 
seen to have a sharp peak at the boundary between Zones I and 
Ila, where it represents a significant fraction of the total heat- 
ing. Outside this narrow heating front, the gas is everywhere 
in approximate thermal equilibrium. The fraction of the bolo- 
metric stellar luminosity that is converted into thermal and 
kinetic energy of the gas can be shown to be — /lad-^oZ-^rad' 
where To (^ 10"^ K) is the Zone I gas temperature and Tiad 
(- Tss) is the color temperature of the incident radiation. For 
the model shown, this fraction is 7%, in good agreement with 
the analytic estimate. 

Other advective models show extremely similar structures. 
The shift in column density of the DF with respect to the static 
models is roughly proportional to A^^, but even models with 
/lad = 0.09 have similar gas temperatures to Model B06 in 
Zones Ila and lib. Models using a Rauch atmosphere instead 
of a black body also produced extremely similar results, de- 
spite this spectrum having a 20 times smaller flux at soft X-ray 
wavelengths (> 54.4 eV). 

A plane geometry is a poor approximation in Zone I for the 
case of photoevaporating knots, in which the ionized flow is 
expected to be transonic and divergent ( jO'Defl et al.||2005| ). 
However, the small observed spatial off'set between the H2, 
Ha, and [N II] emission (OHF07) indicates that the flow in 
Zones II and III is approximately plane parallel. The flow 
timescale through a column density, A^, is equal to N/rtQUo ^ 
32/lad(A^/10'^ cm"^) yr, which is much less than the PN evo- 
lutionary timescale for the H2-emitting portions of the flow, 
justifying our steady-state assumption. However, non-steady 
effects may be important in Zone III, as may the formation of 
CO and magnetic fields, neither of which are included in the 
present models. 

4. PREDICTED H2 SPECTRUM 

Figure |5] shows the radiative efficiency of the models in 
converting the stellar luminosity into H2 emission lines: 77 - 
Liines/if The value of z^tot, corresponding to the sum of all 
H2 lines is < 0.01 for the static models, and rises rapidly 
with /lad for the advective models, approximately as rjiot - 
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Fig. 6. — Excitation diagrams of H2 emission lines for thi'ee Cloudy mod- 
els. Top: Model COO, Middle: Model C06, Bottom: Model A06. Symbols 
show the predicted effective columns from the Cloudy models, determined 
from the emergent line intensities, as a function of the excitation temperature 
of the upper level. Different symbol types indicate the vibrational quantum 
number of the upper level, 02 > as shown in the key. The effective columns are 
normalised with respect to the value for the V2 = Q, J2 = ^ level, which gives 
rise to the 0-0 5 (5) line at 6.907 /jm. Thin lines show the slope of Boltzmann 
distributions for 500 K (solid), 1000 K (dashed), and 2000 K (dotted). Error 
bars show the range of measured values for the Helix knots (see text). 

l-l'lad/Cl + ^^l^) (dashed curve in figure). The dotted line in 
the figure indicates the maximum fluorescent efficiency, as- 
suming that all FUV radiation is converted into H2 lines. 

Figure |6] shows predictions of the full H2 line spec- 
trum for transitions with upper level excitation tempera- 
tures < 17,500 K, which includes most near-infrared and 
mid-infrared lines. Three representative models, are shown 
(Table [TJ: a static model, COO, and two advective models, 
C06 and A06 with, respectively, a low and a high incident 
flux. The line intensities are shown in the standard way as 
effective column densities, which would be proportional to 
exp{-E2/kT) in the case of a Boltzmann distribution at a sin- 
gle temperature, T, giving a straight line on the semi-log plot. 



All static models show a typical fluorescent spectrum with 
strong deviations of the level populations from a Boltzmann 
distribution, whereas advective models show a much smaller 
dispersion in the effective columns of levels with similar ex- 
citation energies, indicating that the excitation is largely ther- 
mal. The slope of the excitation diagram is steeper for lower 
excitation energies, which can be understood with reference 
Figure[T[ the emissivity of the lower pure rotational lines (ex- 
citation temperatures < 5000 K) peaks in Zone lib, where the 
gas temperature is ^ 1000 K, whereas higher excitation lines 
have their peak emissivity in Zone Ila, where the gas temper- 
ature is =^ 2000 K. 

5. DISCUSSION 

Observations of the molecular hydrogen spectrum of the 
knots in the Helix Nebula (Cox e t al | [19981 |H06l |OHF07| l 
are indicated as error bars on Figure |6| It can be seen that 
the two advective models are in broad agreement with the ob- 
servations, whereas the static model is not. Model C06 best 
matches the distance of the spectroscopically observed knots 
from the ionizing star, and indeed shows a better agreement 
than Model A06, which corresponds to the closer-in knots. 
The observed nebular flux in the 1-0 S{1) line and in the sum 
of the 0-0 S(l) to S(7) lines are ^ 1% and - 4%, respec- 
tively, of the bolometric stellar flux ( |OHF07[ l. From Figure |5] 
this can be satisfied by our models with /lad > 0.3//, where / 
is the area covering fraction of knots. A strong prediction of 
our model is that higher excitation lines from the y > 1 levels 
should show a higher population temperature of ^ 2000 K. 
A recent study of an inner knot (Matsuura et al. 2008 ) finds 
a population temperature of about 1800 K for these levels, in 
agreement with the prediction of our relevant model (A06). 

Ro-vibrationally warm H2 has been detected in other PNe 
(e.g., Likkel et al.|2006||McCandhss"eral. 2007) and has fre- 
quently been interpreted as evidence for shocks ( Zuckerman 
& Gatley|1988) . However, EUV-dominated advective PDRs 
may be a promising alternative in these cases too. 
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